Application of hyperbolic scaling for calculation of reaction-subdiffusion front 

propagation 
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A technique of hyperbolic scaling is applied to calculate a reaction front velocity in an irreversible 
autocatalytic conversion reaction A + B — > 2A under subdiffusion. The method, based on the geo- 
metric optics approach is a technically elegant observation of the propagation front failure obtained 
in Phys. Rev. E 78, 011128 (2008). 
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Introduction- The problem of front propagation in 
reaction-transport equations is attracting much attention 
that is related to the considerable progress in our under- 
standing of this phenomenon via the generalization of 
the standard reaction-diffusion scheme in the framework 
of the Fisher-Kolmogorov-Petrovskii-Piskunov (FKPP) 
equation for fractional reaction-subdiffusion systems [l|- 
lll| . The description of reactions under subdiffusion is rel- 
evant to strongly inhomogeneous environments, in porous 
media such as certain geological formations or gels, in 
crowded cell interiors, and so on. Another success in this 
field relates to developing an appropriate new technique 
of treating front propagation, where an appropriate hy- 
perbolic scaling of the reaction-transport equation makes 
it possible to estimate the overall rate of the spreading 
reaction wave without resolving its shape [l2l . [l3j . The 
method of hyperbolic scaling is based on the introduction 
of a small parameter e — > 0, and rescaling of coordinates 
and time (x,t) — > (x/e,t/e), and the contaminant's den- 
sity distribution function. In this case, the problem of the 
wave propagation reduces to the dynamics of the leading 
edge, or the reaction front. Therefore, one analyzes the 
reaction-transport behavior in the leading edge, where, in 
the long-range and long-time limits, the detailed shape 
of the travelling wave is not important. 

In this Report we demonstrate the hyperbolic scaling 
technique to calculate the reaction front velocity in an ir- 
reversible autocatalytic conversion reaction A + B — > 2 A 
under subdiffusion, which, in the case of normal diffusion, 
is described by the FKPP equation [l4|, HU ■ The present 
result is an alternative Hamilton- Jacobi approach via hy- 
perbolic scaling, which is a more elegant presentation of 
the propagation front failure observed in Ref . Q • 

Recent results show that, contrary to normal diffusion, 
the minimal propagation velocity is zero jfj 0, Q , which 
was interpreted as propagation failure (a general discus- 
sion of this issue can be found in [9]). The main fo- 
cus was on situations when subdiffusion can be modelled 
within a continuous-time random walk (CTRW) scheme 
with a waiting-time probability density function (pdf) 
decaying according to the power law, ip(t) ~ t Q , with 
< a < 1 . Analytical and numerical calculations [1(3] 
corroborated this picture, and in the regime of small re- 



action rates, for which the continuous description applies, 
the front velocity was observed to go as . Crossover 
arguments, presented in [ll| . also support this picture. 

Reaction-transport equation - The FKPP equation de- 
scribes a front propagating into the unstable state in bi- 
molecular autocatalytic conversion A + B — > 2A. Ini- 
tially, the whole system consists of particles of type B. 
The introduction of the A-individuals into some bounded 
spatial domain leads to the propagation of a front of A 
into the i?-domain. A general reaction-transport scheme 
that corresponds to the irreversible A +B —> 2 A reaction 
process can be described by the following generalization 
of the FKPP equation 
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Here B(x,t) is a concentration of particles B with the 
initial condition B(x, t — 0) = B — 1, and the condition 
of the mass conservation is 

A(x,t) = 1 - B(x,t) . 

The time kernel M(r) is determined by the waiting time 
pdf in the Laplace domain ip(u) = C[ip(t)] 



M{u) 



uip{u) 
1 - tp(u) 



(2) 



Hyperbolic scaling.- In the sequel, we are concerned 
with the front propagation of particles/individuals of the 
type A. As mentioned above, to analyze the behavior 
of the leading edge, we use the technique of hyperbolic 
scaling developed in pjj], see also 13 1, where the basic 
idea is that, in the long-range and long-time limit, the 
detailed shape of the travelling wave is not important, 
and the problem of wave propagation corresponds to the 
dynamics of the leading edge or the reaction front. We 
follow the details of the analysis presented in Refs. [13, 
[l9| . It is convenient to rewrite Eq. (JXJ) using the variable 
change in the integration with the memory kernel M(t — 
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t') M(t'). Thus Eq. (JTJ for type A individuals reads 
— = a 2 ^ A{M(0(A(x,t-0-l) 

x exp [ - fc / [A(ar, W 

+ fc[l-AO,t)]A(a;,i)- (3) 

After simple manipulation with the second derivative 
over space, it reads 

d A r 

— = a 2 / M(t') A%t-t') -*Vi4(a:,i-0 
Jo ^ 

x / VA(x,t-t")dt" -k(A(x,t-t') - 1) 



x / AA(x,t-i")rft" + fc 2 (A(x,i-i') - 1) 
Jo 

i r vA(x,t-t")dt"] 2 \ 



x exp 



jfc / [A(x,t~t")}dt" dt' 
Jo J 

+ k[l - A(x,t)]A(x,t) , (4) 

where cr 2 = a 2 /2 and we use the t" — > t — t". The 
objective here is to find the rate of the front propagation 
v without resolving the shape of the travelling waves. We 
use hyperbolic scaling for the coordinates x and time t 

x t 
x — » — , t — > - , 

e £ 

and the rescaled density 

A E (x,t)=A( -.- 
\e e 

We write the density A 6 (x,t) in the exponential form 

G E (x,t) 



A e (x, t) = A exp 



(5) 



where the non- negative function G E (x, t) describes the 
asymptotics of the density function and plays a very im- 
portant role in the theory of front propagation. 

At the next step, we rescale x and t variables in Eq. 
@ to obtain 

dA E /"*/ £ r 

e-^ = <t 2 J o M(t'){AA E (x,t-t') 

- kWA E (x,t~t') f WA E (x,t-t")dt" 

Jo 
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x exp 



[A E (x,t-t")]dt' 



dt' 



+ k[l -A e (x,t)]A e (x,t) 



(6) 



We take into account that at finite times in the limit 
e — > the exponent in the Eq. © tends to unity since 
A E (x, t") — > exponentially fast due to Eq. ((5|). Namely, 



lim exp 

E->-0 



[A E {x,t-t")}dt' 



Derivatives of A E (x,t) yield the following expressions 
d t A E = ~(A /e)(d t G E )exp(-G E /e) , 



AA e 



{A /e 2 ){d x G E ) z - (A /s)(d 2 G E ) exp(-G s / £ ) 



We also take into account that the terms of the order 
of (A E (x,t)) 2 in braces tend to zero faster than and 
disappear in the limit e — > 0. Keeping in mind this limit 
and substituting these expressions in Eq. (j6]), one obtains 
for G e = G E (x, t) 
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(7) 



It follows from (J5)) that as long as the function 
G(x,t) = \imG E (x,t) 

E— >0 



(8) 



is positive, the rescaled density A E (x, t) — > as e — > 0. It 
also follows in this limit that d x G(x,t — et') = d x G(x,t) 
and 



lim 

e-s-0 



G E {x,t) G E (x,t-et') 



= [d t G(x,t)}t' . 



Taking into account this limit expression, we obtain the 
following equation for G(x, t) 
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dt 
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In what follows G{x, t) is the action, or Hamilton's prin- 
ciple function, such that 



H 



dG 
"dt' 



P 



dG 
dx 



(10) 



are the Hamiltonian, and the momentum, respectively. 
Therefore, it follows from Eq. ([TUf that Eq. © is a kind 
of Hamilton- Jacobi equation. The Laplace transform of 
the memory kernel yields 



M{H) = CM(t) / M{t')e 
Jo 



- m 'dt' 



1 - JW) 



(11) 
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where ip(H) = Cip(t) is the Laplace image of the wait- 
ing time pdf [see Eq. (J2J]. Finally, one obtains that 
the Hamiltonian can be found from the Hamilton- Jacobi 
equation 



dG 
~dt 



\~dx~J 



-**(^YM(H) 



kg 2 (dGy 
Eventually, it reads 



M(H) - M(0) 



- k. 



(12) 



H = a 2 p 2 M(H) 



ka^ 

IT 



v 



M(H)-M(0) +k, (13) 



and the action is G(x, t) — f \p(s)x(s) — H(p(s),x(s))]ds. 

The rate v at which the front moves is determined from 
Eq. (J5J) at the condition G(x,t) = 0. Together with the 
Hamilton equations, this yields 



dH 
dp 



H 

P 



(14) 



Note that the first equation reflects the dispersion con- 
dition, while the second one is a result of the asymptoti- 
cally free particle dynamics, when the action is G(x, t) — 
px — Ht. Taking into account x = vt, one obtains Eq. 
(fl4l) (s ee also details of this discussion e.g. in Refs. 
[13L Yw[). Now we analyze these two Eqs. (fT4|) to de- 
fine v. 

Markovian case- 
First let us check the Markovian case, when ip = 
- exp(— t/r). Thus one has 



1>(H) = 



1 



1 + Hr 



where t is a characteristic time scale. Therefore, from 
Eq. (jTTJ) the Laplace image of the memory kernel reads 



M(H) = M(0) = 1/t. 



(15) 



In this case, the Hamiltonian in Eq. (Tl"3"|) is H = Dp 2 + k, 
where a 2 /r = D is a diffusion coefficient. The moment 
is 



p(H) 



(H-k). 



From Eqs. (fT4|) we have 2Dp = H/p and H = 2k, which 
yields for the overall velocity of the front propagation 

v = 2VkD, (16) 

which is the classical FKPP propagation speed (see dis- 
cussions in Refs. d, HH). 

Subdiffusion.- We have for subdiffusion i/j = 1+( ^ t ^i +a , 
which yields @, QJ] 



M(H) 



H 1 



, (M(0) = 0), 



(17) 



where the transport exponent is defined in the range < 
a < 1. The Hamiltonian is 



H = D a (H- k)H- a p 2 



(18) 



where D a = a 2 /r a is a generalized diffusion coefficient 
From Eq. (fT8|) one obtains 




and Eqs. (fT4")) result in 



op p a 



(19) 



(20) 



This equation has the solution for H = H = 0. There- 
fore, the asymptotic velocity of the front propagation is 

v(H) = 0. 

Conclusion.- We demonstrated the technique of hyper- 
bolic scaling for the calculation of the reaction front ve- 
locity in an irreversible autocatalytic conversion reaction 
A + B — > 2A under subdiffusion. This is a technical 
presentation of the powerful Hamilton- Jacobi method for 
asymptotic estimatimation of the propagation front ve- 
locity observed in Ref. @ ■ 

It should be admitted that the dispersion velocity v(H) 
also determines the relaxation rate at the large time 
asymptotic t — > 00 for the finite value of H. A qualitative 
crossover argument based on the truncated power-law 
distribution was suggested in Ref. [llj . According to the 
arguments, for short times the behavior of the velocity 
v (H) must be similar to that in subdiffusion (it does not 
feel the cutoff), whereas for long times the behavior is the 
classical one with a constant minimal velocity, and there 
has to be a crossover (no jump!) at a time t CI between 
the two of them. Assuming that the time dependence of 
the velocity in the anomalous domain is v (t) ~ t and, 
after the crossover to the normal domain, the velocity, 
determined by Eq. (|16l) . v = 2\JkD is attained, both can 
be equated at t cr to obtain the corresponding j3. To de- 
termine the crossover time, it is plausible to argue with 
the amount of performed steps, as a measure of mobility, 
which for the normal regime is nu{i) cx t/t 1 ^ a T a and in 
the subdiffusive regime reads nso(t) oc (t/to) a , equating 
them, one finds v(t < t a ) ~ . 

The hyperbolic scaling also corroborates the relaxation 
picture for the velocity v(H) ~ i~ obtained in Refs. 
[lOl fll| . An important point when considering the re- 
laxation in the framework of hyperbolic scaling is that 
the process of relaxation for a subdiffusive front can be 
treated for the finite energy H = 2k in the framework 
of the Markovian case. For normal diffusion, the hyper- 
bolic scaling method is rigourously justified [12J, and the 
method yields a correct result for v(H) in Eq. (I16|) . which 
is exactly the FKPP case. This can be demonstrated for 
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a truncated waiting time pdf ipT{t)- The latter is conve- 
nient to take, e.g., in the following power-law form 



1 + (t /T) a 



,-t/T 



i + (t/t r+ i 



(21) 



where T has the role of the cutoff. Therefore, for any 
finite T the mean waiting time is finite: 



at? 



1 + (to/T)« 



T 



For this normal diffusion hyperbolic scaling yields Eq. 
(jT6|l for the velocity in the form 



v(T) oc VDkxT^- , 



(22) 



which corresponds to the relaxation rate obtained in 
Ref. [ilj . 

Another specific property of the method is an effective 
linearization of the generalized FKPP Eq. (j4} . It should 
be admitted that this relates to considering a wavefront 



with an exponentially decaying leading edge moving with 
a constant velocity v. The standard, traditional way to 
perform this analysis is first to linearize the equations. 
Hyperbolic scaling performs it automatically, since the 
density A is not zero only when G = 0. Moreover, it 
also affects the integrand kernel in Eq. (j4]), namely, as 
admitted above, in the limit e — > 0, the exponent in the 
Eq. © tends to unity since A 6 (x,t") — > exponen- 
tially fast due to Eq. ([5]). This essential simplification 
makes it possible to apply the strong machinery of the 
Laplace transform and arrive at the analytically treat- 
able Hamilton- Jacobi equation (fT2"j) that is an easy and 
elegant way to obtain the front propagation, namely the 
failure of the latter. This nonlinear kernel was also stud- 
ied in relation to a mechanism coupling the waiting time 
distributions to the reaction [7j to resolve a controversy 
about reaction-subdiffusion front propagation. To this 
end, a more general scheme of the local waiting time 
was suggested [7j that eventually leads to a more compli- 
cated analysis in the framework of the Hamilton- Jacobi 
approach than presented here. 
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